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Abstract. We use semi-analytic modelling of the galaxy-cluster population and its strong lensing efficiency to explore how the 
expected abundance of large gravitational arcs on the sky depends on cr 8 . Our models take all effects into account that have been 
shown to affect strong cluster lensing substantially, in particular cluster asymmetry, substructure, merging, and variations in the 
central density concentrations. We show that the optical depth for long and thin arcs increases by approximately one order of 
magnitude when cr 8 increases from 0.7 to 0.9, owing to a constructive combination of several effects. Models with high cr 8 are 
also several orders of magnitude more efficient in producing arcs at intermediate and high redshifts. Finally, we use realistic 
source number counts to quantitatively predict the total number of arcs brighter than several magnitude limits in the R and I 
bands. We confirm that, while <x 8 ~ 0.9 may come close to the known abundance of arcs, even <x 8 ~ 0.8 falls short by almost an 
order of magnitude in reproducing known counts. We conclude that, should cr 8 ~ 0.8 be confirmed, we would fail to understand 
the strong-lensing efficiency of the galaxy cluster population, and in particular the abundance of arcs in high-redshift clusters. 
We argue that early-dark energy or non-Gaussian density fluctuations may indicate one way out of this problem. 
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1. Introduction 



After the 3-year WMAP data release dSpergel et alj|2007l) . the 
best-fit value for the normalisation of the density-fluctuation 
power spectrum was lowered from cr 8 ~ 0.9 to cr 8 ~ 0.75. This 
puts many of the cosmological tests based on structure forma- 
tion under stress, su ch as the number counts of galaxy clus- 
ters and its evolution dEvrard et al.l2008 ; Rines et alj2007[) . the 
large-scale structure probed via weak gravitational lensing, and 
large optical surveys that t end to favour a value of o~% ~ 0.9 as 
well (Hoekstra et al .I2006I) . The value of erg is perhaps the least 
well-known of the main cosmological parameters today. The 
controversy is highlighted by the fact that largely discrepant 



value s of cr% are being published dReiprichll2006t lEvrard etal 
2008) 



At the same time, there are signs of convergence on a nor- 
malisation which may be compatible with most measurements 
of the amplitude of structures large enough for linear evolution 
to dominate. At <x 8 ~ 0.8, CMB and gravitational-lensing mea- 
surements may meet after a slight increase in the mean redshift 
of the background so urces used to identify the weak-lensing 
signal JFu et alJ2 008). Some analyses of the X -ray cluster pop - 
ulation also seem to find this value acceptable (Reiprich 2006). 

Suppose, then, that <x g ~ 0.8. We argue here that this ex- 
acerbates a problem with non-linear structure growth on the 
cluster scale to a level which seems serious despite consider- 
able uncertainties. Specifically, we shall address the statistics 
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of gravitational arcs in galaxy clusters to show that none of the 
many possible explanations suggested in the past decade suf- 
fices to bring theoretical expectations into agreement even with 
the admittedly sparse observations. 

We summarise the situation in Sect. 2 and compile some 
analytical results concerning model universes with different 
normalisation of the power spectrum in Sect. 3, in order 
to gain quantitative insight into the situation. In Sect. 4, 
we describe the semi-analytic modelling of lensing by the 
galaxy-cluster population. Section 4 summarises the results, 
and Sect. 5 presents our conclusions. There, we also discuss 
how our results comply with the new WM AP-5 data release 
( Dunklev et al.l I2008I : iKomatsu et all I2008I) which was pub- 
lished after this work was completed. 

2. Is there an arc-statistics problem? 

Comparing the observed abundance of gravitational arcs to the 
efficiency of numericall y simulated galaxy cluste rs for strongly 
lensing distant galaxies, iBartelmann et al. I dl998l) claimed that 
about an order of magnitude fewer arcs are expected in the 
ACDM cosmology than are actually observed. This was called 
the arc-statistics problem. Numerous attempts were carried out 
to see how the problem could be solved. A corrugation of the 
lensing potential by individual cluster galaxies turned out to 
have a slight net effect on large arcs because the increased 
length of the causti c curves was counteracted by splitting arc s 
into shorter pieces (lFlores et al. 2000; Meneghetti et alj|2000l) . 
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Cluster asymmetry, however, was ident ified as crucial 



for the abundant formation of arcs (Bartelman n et al.l 11995 



Molikawa & Hattorill200U iMeneghetfi et al.ll2007l) . Simplified 
analytic clu ster models were unable to reproduce the nu mer- 
ical results (ICoorav Il999t Kaufmann & Straumann l2000l) for 
two reasons; first, they did not account for the dependence of 
cluster concentrations on the cosmological constant; and sec- 
ond, simple elliptical cluster models were also found inad- 
equate to quantitatively explain th e arc-formation efficienc y 



in numerically simulated models (Menegh etti et al.l 12003b) 



Lensing properties of numerically simulated galaxy clusters 
were found to agree we ll with those of comparable, real clus- 
ters dHoresh et al.ll2005h . 



Based on numerical simulations. IWambsganss et al.l (12004) 
found that the probability for high magnifications along light 
rays propagating to us from sources in the distant Universe de- 
pends steeply on the source redshift distribution. Based on this 
result, they concluded that there is in fact no arc-statistics prob- 
lem if the realistically distant tail of the source distribution is 
taken into account. However, their identification of highly mag- 
ni fied light bundle s with strongly distorted arcs was questioned 
bv lLi et al.l (12005b . who showed that at least numerically simu- 
lated clusters produce a class of highly magnified, but weakly 
distorted images. In particular, numerical simulations show that 
the most likely length-to-width ratio for arcs magnified by a 
factor > 10 is typically ~ 3 (see also Figure [9] below). High 
magnification probability does therefore not imply a frequent 
occurrence of strongly distorted arcs. Identi fying arcs in sim 



ulations, the redshift dependence found by IWambsganss et al 



(2004) was qualitatively confirmed, but quantitatively fou nd to 
be substantially weaker dLi et alfeOoltlFedeli et al.ll2006h . 
The arc-statistics problem was again questioned by 



Dalai et a l. (2004), who used nu merical simulations to con firm 
the overall optical depth found by Bartelmann et al.l 0998), but 
found agreement between the observed and expected numbers 
of arcs because they estimated a lower abundance of observed 
arcs and a higher number density of background sources. They 
took the sou rce-redshift distribution in to account, but found it 
weaker than Wambsganss et al. (2004) had claimed. 

The problem did not disappear, however, because it 
was observationally found that in particular the number 
of arcs in distant clusters is considerably higher than 



naively expect ed (Thompson et a l.11200 It Gladders et al.1 12003. 
Zaritskv & Gonzalezll2003l) . IOguri et aLT(l2003b studied the ef- 



fect of central concentration and triaxiality of cluster-sized ha- 
los on their strong-lensing ability, finding that triaxial clusters 
with sufficiently steep central density profile can come close to 
the observed results. While this effect is undoubtedly present, 
it is also included in numerically simulated clusters and does 
thus not attenuate the apparent discrepancy between the arc 
abundances observed in the sky and produced in realistically 
simulated cluster populations. 

So, is there an arc-statistics problem? Based on the preced- 
ing discussion, the question still seems undecided. However, it 
is important to note that up to this point, all simu lations and cal- 
culations were done assuming 0.9 < cr% < 1.2. lLi et al. (2006) 



analysis, cr 8 = 0.74 dSpergel et al.l2007l) . lLi etai] d2006l) found 
that the expected abundance of large arcs on the sky drops very 
steeply if erg is lowered, reflecting the exponential decrease of 
the abundance of massive clusters with aj,. 

Motivated by this earlier study, in which numerical simula- 
tions with two different values of <x 8 were compared, we shall 
investigate in this paper how the number of gravitational arcs 
predicted in a standard cosmological model changes with crs, 
and how the results compare with the observed statistics. Such 
a study is possible only because we can replace numerical sim- 
ulations by se mi-analytic calcula tions using a novel algorithm 
developed bv lFedeli et al. (2006). The method captures two in- 
gredients crucially important for stro ng cluster lensin g: clus- 
ter asymmetries and cluster mergers dTorri et al.ll2004l) . which 
were fou nd to have transien t, but strong effects on arc cross 
sections ( Fede li et al . 2006). We also properly account for a 
realistic source redshift distribution and the luminosity func- 
tion of ba ckground source s . The c ontribution of cD galaxies is 
ignored. Meneghet ti et al. ( 2003ah showed that the cross sec- 
tion for giant arcs of numerical clusters is increased by < 50% 
when a realistic cD model is included, and this does not signifi- 
cantly alter our conclusions. In addition, the same paper shows 
that the boosting effect of a massive central galaxy is much 
reduced when the host model cluster is asymmetric, as those 
used here. We also neglect the effects o f gas physics, which are 
potentially important in cluster cores. IPuchwein et al.l (|2005) 
have shown that, depending on the physical effects included 
into the modelling, the lensing cross section of individual clus- 
ters can be i ncreased by < 100% due to the presence of bary- 
onic matter. IWambsganss et al.l d2008) implemented a specific 
model describing baryon cooling and galaxy formation in the 
core of cluster-sized dark matter halos, finding an in crease in 
the production of a rcs of ~ 25%. Hilbert et al.l (12008) (see also 
Hilb ert et al.l20 07) performed a very similar study on the effect 



pointed out that the problem is substantially aggravated if erg 
is lowered to the value preferred by the third-year WMAP data 



of the stellar component on optical depths for image splitting 
in the Millennium Simu lation, finding re sults compatible with 
earlier analyses dMeneghetti et al.1 l2000h . Although baryonic 
physics may increase strong-lensing cross sections, it is not ex- 
pected to bridge the order-of-magnitude gap between theory 
and expectations, although the uncertainties in the modelling 
are admittedly substantial. 

We shall consider five different cosmological models. The 
energy density content and the Hubble constant are taken 
from the WMAP- 3 data combined with the SDSS observations 
dSpergel et al. Il2007l) . and kept fixed throughout the paper. They 
are Q m o = 0.265, Qa,o = 0.735 and h = 0.71. The values of 
erg are chosen differently for the five models and are 0.7, 0.75, 
0.8, 0.85 and 0.9 respectively. This choice allows us to cover 
the complete range of values from the third-year WMAP data 
to normalisation values derived from galaxy-cluster counts. 
Typical degeneracies among different cosmological parameters 
imply that changes in <r% would cause O m and possibly other 
parameters to change, depending on the data set underlying the 
parameter determination. However, since the main purpose of 
this work is to isolate the impact of the power-spectrum nor- 
malization on arc statistics, we choose to keep everything fixed 
except erg. We shall discuss the effect of relaxing this assump- 
tion in the last section of the paper. 
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Fig. 1. Left pane l: Conc entration of dark matter halos accord- 



ing to lEke et alJ Q200 1 ) as a function of virial mass at fixed 



redshift, z - 0.3. Five different values for <t& are considered, 
as lab elled in the plot. R ight panel: Collapse redshift accord- 



ing to lEke et al.l yOOl) as a function of <x 8 for dark-matter 
halos at z = with three different masses, M = 10 14 M o /i _1 
(black solid line), M = 5 x l0 u M o h~ 1 (blue dashed line) and 
M = IO^Mq/T 1 (red dot-dashed line). 
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Fig. 2. Left panel. iPress & Schechterl d 1 974 ) mass function at 
fixed redshift z — 0.3 as a function of mass. Right panel. Merger 
rate between a cluster of mass M - 10 15 M Q h~ l and a substruc- 
ture of mass m indicated on the abscissa, at z — 0.3. Both pan- 
els show results for different values of <xg, as labelled. Note 
that here and in Figure [3] the merger rate is the probability for 
a dark-matter halo to merge with a substructure per unit loga- 
rithm of the merging mass and per unit logarithmic time. 



3. Expectations 

We need theoretical predictions of the arc abundance which are 
at the same time as precise as possible and fast to achieve be- 
cause they need to be carried out in many cosmological models. 
We first try to gain some insight into the various contributions 
to the arc optical depth, its dependence on and its variation with 
the normalisation of the power spectrum. 

We start with the internal structure of dark matter ha- 



los that, as described in a variety of s t udies ( Cole & Lacev 



19961: INavarro et all Il997t llmgl boOOt iBullock et all 12001 
Eke et al-lkOoilT depends on the complete halo formation his 



tory. Assuming that the den sity profile of cluster-sized dar k 
matter halos is of NFW form dNavarro et al.ll995Ul996lll997h . 
we first explore how the concentration of the profile depends on 
cr 8 . The concentration of a dark matter halo is the ratio of the 
virial radius to the scale radius r s of the density profile. Here 
we define the virial radius as the radius of the sphere inside 
which the average density of the halo is 200 times the critical 
density of the Universe at the given redshift. This radius sepa- 
rates well the internal, relaxed part of isolated galaxy clusters in 



numerical simulations from the external, infall part (lEke et al 



1998). We expect that higher <r 8 allows earlier structure forma- 
tion, such that halos form from a higher background density 
and have more time to relax, causing higher concentrations. 

This expectation is verified in the left panel of Figure Q] 
w here the halo con centration as obtained from the prescription 
bv lEke etaL ( 200 lb is shown as a function of halo mass at fixed 



redshift z = 0.3, typical for strongly lensing clusters. Apart 
from the well-known decrease of concentration with halo mass, 
the figure shows that the concentration tends to increase with 
the normalisation. For example, the concentration of a cluster- 
sized dark-matter halo of M — 1O 15 M /i _1 increases by ~ 40% 
when erg is increased from 0.7 to 0.9. This fact alone may sig- 
nificantly affect arc statistics, because more compact cluster 
cores push the critical curves and caustics outwards, thus in- 
creasing their strong-lensing cross sections. 



According to lEke et al.l (120011) . the collapse redshift z c of a 
halo of mass M is implicitly given by 



D + (z c )(t(M s ) 



d In cr(M s ) 



d\nM 



1 

C ' 



(1) 



where D+ is the linear growth factor for density fluctuations, 
normalised to unity at present, and M s is the mass contained 
within the radius of maximum circular velocity for the NFW 
density profile, 2.17r s . The rms density fluctuation <j(M) is 
taken at the linear scale corresponding to M, and C = 28 is a 
di mensionless cons tant calibrated against iV-body simulations 
by lEkeet allfeOOlh . 

Evidently, higher <x 8 implies a lower growth factor at col- 
lapse redshift in equation (HJ, hence a higher collapse redshift. 
We show the collapse redshift for dark matter halos of different 
mass at redshift zero as a function of the cr 8 in the right panel 
of Figure Q] 

Alternative prescriptions for the computation of halo con- 
centrations exist, but affect arc statistics onl y mildly. For 



instance, 



Fedeli et all d2007l) s howed that the INavarro et al 



(1997) and Bulloc k et al 



2001) recipes overestimate and un- 
derestimate respectively the cro ss section for g iant arcs by a 
factor of ~ 2 with respect to the lEke et al.l (1200 11) prescription. 
However the latter is probably the most general and physically 
best motivated. It turned out to reproduce halo concentrations 
in a variety of cosmologies, i ncluding those with dynamical 



in a variety or cosmologies, in 
dark energy (foolag et al.ll20 04). 

Next, we analyse the halo mas s function and the merger 



rate. Figures [2] and [3] show the Press _& Schechter (1974) 



function and the merger rate dLacey & Cole 



1993 



1 19941) at fixed redshift as a function of mass, and at fixed 
mass as a functi on of redshift, respectively. We neglect 



the improvements ( Jen kins et al.ll2001t ISheth & Torm en 2002 



Warren et al.l l2006) of the mass function here because we only 
mean to illustrate the differences between different normali- 
sations. The behaviour of the mass function is quite obvious. 
Higher normalisation gives rise to more structures at a given 
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Fig. 3. Left panel. Halo number density obtained from the 
Press & Schechterl (119741) mass function as a function of red- 
shift at a fixed halo mass M = 7.5 x 10 u M Q h- 1 . Right panel. 
Merger rate between a cluster of mass M = 10 l5 M o h~ l and a 
substructure of mass m = 5 X 10 13 M e h~ l as a function of red- 
shift. Again, results are shown for five different values of <x 8 . 



redshift. This is particularly evident at the high-mass end where 
the mass function depends exponentially on cr^. At the massive 
cluster scale, the mass function can vary by orders of magni- 
tude as cr 8 is varied. 

The behaviour of the merger rate may be less obvious. It 
is larger for higher erg if the mass of the main halo is large 
and the mass of the secondary halo is a considerable fraction 
of it. These are rare events because massive structures (and 
substructures) are rare. On the other hand, substructures much 
less massive than the main halo merge at a lower rate in highly 
normalised models. Because of the much higher abundance of 
low-mass halos, we expect mergers with low-mass substruc- 
tures to be equally or even more frequent in models with lower 
<t 8 , at least in the redshift interval relevant for our purposes. 
Since strong cluster lensing is highly sensitive to asymmetries 
and external perturbations, more merger events will further in- 
crease the strong-lensing cross sections. 

Observed strong-lensing clusters frequently show substruc- 
tures. While this could introduce a bias on the observational 
side, the study we s hall refer to in this w ork is based on X- 
ray selected clusters ( Le Fevre et al.|[l994l) . An investigation of 
the possible selection effects introduced by this fact, account- 
ing for the boosting effect that cluster mergers have on both 
the lensing efficiency and the X-ray luminosity, is reported in 
Fedeli & Bartelmannl d2007bl) . We are preparing a study of the 
strong-lensing properties of merging clusters in a large cosmo- 
logical simulation. 

The larger number-density of halos, their higher concentra- 
tions and the modified merger activity with increasing <x 8 will 
all combine to make the strong-lensing optical depth depend 
sensitively on the normalisation. It will be critically important 
for precise predictions to include the effect of mergers into the 
calculation. We shall now present our results, confirming these 
expectations. 



4. Cluster population 

To produce realistic models of the cluster population with- 
out time-consuming numerical simulations, we adopt the ex- 



tended Press-Sch e chter formalis m ( Press & Schechterl 1 1974c 
Bond et al. 1991 : Lacev & Cole 1993b to construct merger 
and formation histories for a set of Af = 1 , 000 dark-matter 
halos for each model universe. E xamples for t he applica- 
tion of this procedure are g iven in iRandall et ; 
Fedeli & Bartelmannl d2007al) . among others. 
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To cover the mass range relevant for strong cluster lens- 
ing, we draw the dark-matter halos uniformely from the mass 
interval between 10 14 M o /z _1 and 2.5 x IO^Mq/T 1 at redshift 
z — 0. Each object is then evolved backwards in time in suitably 
chosen discrete time steps. Lensing by each dark-matter halo 
at each redshift step is modelled using a NFW density profile 
whose lensing potential is elliptically deformed with an ellip- 
ticity of e = 0.3. This value was shown to give the best agree- 
ment with deflection ang le maps of realistic simulated clusters 
dMeneghetti et alj|2003bi) . 

The transient boost of the strong-lensing efficiency due to 
cluster mergers is taken into account by modelling the merg- 
ing substructures also as elliptical NFW lenses. Each time a 
main cluster halo undergoes a merger with a substructure hav- 
ing more than 5% of its mass, the encounter is modelled assum- 
ing that the two halos approach at a constant velocity starting 
from a distance equal to the sum of their virial radii, and as- 
suming that the merger proceeds at the gravitational free-fall 
time. 

Given the deflection-angle maps for each model cluster at 
each redshift, we compute the strong-lensing cross section for 
arcs with length-to- width ratio d > 7.5 and d > 10. Source 
redshifts were drawn randomly from the distribution 



P(Zs) = 



P 



z; exp 



(2) 



dSmail et al.ll 19951) . with z = 1 and /3 = 3/2. The distribution 
peaks at z ~ 1.2, implying that objects at z — 0.3 - 0.5 are the 
most efficient lenses. Model clusters were evolved backwards 
in time up to the source redshift. 

Phot ometric redshif t s measured for recent wide-area sur- 
veys dllbert et al. 2006; Se mboloni et al. 20061) favour some- 
what different values for the parameters of equation © or even 
a different functional form for the redshift distribution, imply- 
ing a peak at z < 0.7. While such distributions would be signifi- 
cantly less efficient in producing large amounts of long and thin 
arcs, it is likely that these wide surveys are not deep enough to 
capture the complete source population relevant for the produc- 
tion of strong-lensing features. This is suggested by the high 
mea sured redshifts of man y giant arcs in galaxy clusters (see 
e.g. lElfasdottir et al]|2007l for a recent example). 

We compute the lensing c ross section s by m eans of the 
semi-analytic prescription by iFedeli et al.l (120061) instead of 
costly ray-tracing simulations. The method consists of integrat- 
ing the inverse magnification along critical curves over an area 
determined by the locally linearised lens mapping. In particu- 
lar, the area is defined by the condition \A I /A t \ > d, where A r 
and A t are the (radial and tangential) eigenvalues of the lensing 
Jacobian matrix. Finite source size is accounted for by con- 
volving the lens properties with functions standing for circular 
background galaxies with 0.5" radius, and source ellipticity is 
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Fig. 4. The optical depth per unit redshift for arcs with length-to-width ratio d > 7.5 for five different values of cr 8 , as labelled 
in the plots. The left and right panels show, on the same scale, results including and excluding the effect of cluster mergers, 
respectively. Predictions for d > 10 are extremely similar and therefore not shown here. 



in cluded in the co mputation by using the f ormalism develope d 
by lKeetonl (120011) . Further detail is given in lFedeli et al.l (120061) . 

Having computed the strong-lensing cross section of each 
individual cluster at each redshift step, we can compute the dif- 
ferential optical depth per unit redshift dfj/dz. As mentioned, 
we use the common thresholds d = 7.5 and d — 10 for the 
length-to-width ratio, but show results for only one of them if 
the predictions have similar behavior. Since the cluster popula- 
tion is represented by a discrete set of objects, 

dfafe) y 1 g-d(M,-,z,z M ) f M ' +1 .... ... 

—a — = / , a — ^2 n(M,z)dM , (3) 

where D s , is the angular-diameter distance to the source sphere 
of the r'-th cluster, and n(M, z) is the number of structures with 
mass within M and M+dM, and redshift within z and z+dz. The 
masses M, are assumed to be in ascending order, M;+i > M, 
for all i. Since source redshifts are assigned randomly to each 
cluster in our synthetic sample, the weighting with the redshift 
distribution is implicitely included in the calculation. 

The integral over the lens redshift of equation (0 is the av- 
erage optical depth. After multiplication with the total number 
of available sources in the sky, it yields the total number of 
gravitational arcs with length-to- width ratio > d, 

N d = n s f d = n s P ^dz . (4) 
Jo dz 

It should be noted that the integral on the r.h.s. of equation (|4]i 
cannot extend to infinity for our discrete cluster sample, but 
only to a finite value z max where the probability of finding 
sources according to the adopted redshift distribution equa- 
tion (f2]i is negligib le. We adopted Zmax = 7.5 (see also the dis- 
cussion in Fedeli & Bartelmann 2007a). 
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Fig. 5. The total average optical depth for arcs with length-to- 
width ratio d > 7.5 as a function of erg, The red dot-dashed line 
includes mergers, while the black solid line ignores them. 

5. Results 

5.1. Optical Depth 

Figure [4] illustrates the strong-lensing effects of our synthetic 
cluster population. It shows the optical depth per unit redshift 
given by equation (O, computed for the five different normal- 
isations <t 8 analysed here. We also show the change in optical 
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Fig. 6. The ratio of the total average optical depth for arcs with 
length to width ratio larger than or equal to d — 7.5 (black solid 
line) and d — 10 (red dot-dashed line) to that obtained ignoring 
the boosting effect of cluster mergers, as a function of erg. 

depths caused by cluster mergers. Lines are drawn only for one 
of our choices of the length-to-width threshold for gravitational 
arcs, namely d = 7.5. Differential optical depths for d = 10 are 
slightly smaller, but the qualitative behavior is unchanged. 

As anticipated in Sect. [3] the optical depth per unit redshift 
increases substantially with <x 8 because more halos with more 
concentrated mass distributions exist as <x 8 increases. This dif- 
ference is particularly strong at high redshifts, z > 1, where the 
lensing efficiency is still significant for <x 8 = 0.9 while being 
negligible for <x 8 = 0.7. In this sense, a higher normalisation is 
degenerate with the introduc tion of an early-dark energy com- 



T 

ponent (see the discussion in Fedeli & Bartelmann 2007a b 



Figure [5] shows the total optical depth according to equa- 
tion dUi, i.e. the integral over each of the curves shown in 
Figure H] for five values of cr 8 , either including or neglecting 
cluster mergers. There is more than one order of magnitude dif- 
ference between the model with lowest and the highest normal- 
isations, <x 8 = 0.7 and <x 8 = 0.9, independent of whether clus- 
ter mergers are included. This high sensitivity of the expected 
number of arcs on erg has important consequences for the arc 
statistics problem, as will be discussed in detail in Sect. [6] 

Figure |5]also shows that mergers are more important when 
erg is low. This is emphasised in Figure [6] which shows the 
ratio between the total optical depth obtained including and 
neglecting cluster mergers. Evidently, dynamical activity en- 
hances the optical depth by a factor > 5 if <x 8 is low, and by 
a factor > 3 when cr 8 is high. This co nfirms results obtained 
earlier bv iFedeli & Bartelmann! d2007al) . We also note that the 
effect of mergers is high for more extreme arcs, d > 10. This 
is because the individual cross sections are smaller in this case 
and therefore relatively more sensitive to perturbations. 



It is also highly interesting to see how the strong-lensing 
efficiency of clusters at high redshift changes with cr 8 . This 
is particularly relevant in view of the apparent unexpectedly 
high incidence of gravitational arcs in distant g alaxy clusters 
(bladders et al.l2003tlza"ritsky & GonzaleJ2003l) . We quantify 
this by means of the cumulative optical depth 



C d (z) 



■r 



df d (z') 
dz' 



dz' 



(5) 



contributed by clusters above redshift z (the upper integration 
limit is set to z max = 7.5 here as well). Figure|7]shows Cd(z) for 
different normalisations <x 8 relative to Cd(z) for <x 8 = 0.7, in- 
cluding and neglecting mergers. This ratio is virtually indepen- 
dent of the length-to-width threshold, so we show it for d > 10 
only. 

The relative contribution to the optical depth increases to- 
wards high redshifts for all models with <x 8 > 0.7. The increase 
is higher when mergers are ignored (~ 3 instead of ~ 2 orders 
of magnitude) because the relative effect of mergers decreases 
as unperturbed clusters by themselves become stronger lenses. 
The main conclusion is that high <x 8 makes gravitational arcs 
at high redshift extremely more likely. Conversely, this means 
that arcs in distant clusters are a massive problem for models 
with low <x 8 . 

5.2. Number of Arcs 

We now need to transform our predictions for the optical depth, 
including mergers and a realistic source redshift distribution, 
into predictions for the numbers of observable arcs on the sky. 
To do this, it suffices in principle to multiply the average optical 
depth f d with the total number of sources in the sky, according 
to equation (|4). 

However, it is necessary to include a luminosity function 
for the faint blue background galaxy population into the calcu- 
lation and to account for the magnification effect due to grav- 
itational lensing. The latter has a twofold impact on the ob- 
served source counts. First, faint sources are magnified above 
the flux threshold for detection, thus increasing the number 
of sources visible per unit solid angle. Second, the sky is lo- 
cally stretched, thus reducing the number density of the source 
galaxies. If the original flux distribution function of the sources 
is a power law with logarithmic slope -1, then the two effects 
of flux magnification and number dilution cancel, leaving the 
number of observed sources p er unit solid angle unchanged 
dBartelmann & SchneideifeoOll) . 

We model the number counts of faint background galaxies 
as a function of t he observed appa rent magnitude to match the 



measurements of iCasertano et al.l (120001) in the Hubble Deep 



Field. There, only the I-band magnitudes are used, and the 
number counts are fitted using the relation 



n (mi) = 



«o,* 



sJlQlaimu-ii'i) + 102i(mi,i-mi) 



(6) 



with the best-fit parameters a = 0.30, b = 0.56, m^i = 20 
and no,* = 3 x 10 3 per square degree and per unit magnitude. 
Since these number counts are given in the I-band only, we 
convert them to the R-band by using the approximate relation 
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Fig. 7. The cumulative optical depth for different cr% is shown relative to its value for <x 8 = 0.7. The left and right panels show 
results including and neglecting cluster mergers, respectively. The results are quite independent of the length-to-width threshold; 
d > 10 was chosen here. 



R ^ I + 1. In Figure |8] the colour R - I is shown as a function 
of redshift for different morphological types of galaxy. Since 
gala xies imaged as long and thin arcs are usually blue spirals 
(see lTyson & Seitzerll98ilEliisl 19971) . R-I ~ 1 approximately 
holds over the redshift range z < 1.7. However, also adopting 
the relation R — I ~ 0.5, which holds better at higher redshift, 
the magnified number counts (see below) change by less than 
50%, which is irrelevant given the un certainties in the original 
galaxy counts (ICasertano et al.1 1200 0). We thus conclude that 
the approximation R I + 1 is admittedly rough, but sufficient 
for our purposes. 

Figure [9] shows the conditional probability distribution for 
the magnification of background sources given the length-to- 
width threshold for the imaged arcs. The results shown are ob- 
tained from fully numerical ray-tracing simulations, to which 
we fit two-component Gaussians for both thresholds d — 7.5 
and d — 10. The two-component Gaussian is 



p(mM = 



V27TO-! 

1 -A 



exp 



2cr\ 



V27TI 



exp 



7TCT2 



2cr 2 2 



(7) 



with ju + = \fi\. The best-fit parameters A, <j; and f/.+j are sum- 
marised in Tab. Q]for both d = 7.5 and d = 10. Obviously, the 
magnification is not a good estimator for the length-to width 
ratio of an image, especially when high length-to- width thresh- 
olds are used. 

The original number counts read off Casertano et al. 
are thus convolved with the conditional probability distribution 
for arcs with thresholds d = 7.5 and d — 10 to obtain the num- 
ber counts after the magnification bias. Let n^(F) be the origi- 
nal flux distribution function for the sources, that is the number 
of sources per unit solid angle contained in the unit flux inter- 



2.5 



h — i — i — i — I — i — i — r 




Fig. 8. The R-I colour as a function of redshift, computed from 
the spectra of three different morphological types of spirals and 
for elliptical galaxies, as labelled in the plot. 

val around F. It corresponds to the magnitude distribution of 
equation (O with the magnitude replaced by the flux. Then the 
magnified distribution is 



f +M (F\ P( f i + \d) J 

n(F) = « — — dfi + . 

Jo \f* + 1 H% 



(8) 



The magnified number counts can then simply be multiplied 
with the average optical depth to find the total number of 
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Fig. 9. Conditional probability distribution for the magnifica- 
tion of images given a threshold for the length-to-width ratio 
of d — 7.5 (left panel) and d — 10 (right panel). The black 
solid lines show the result of ray-tracing simulations, while 
the red dot-dashed curves are the best-fitting two-component 
Gaussians whose parameters are summarised in Tab. Q] 



arcs. Figure \TU\ shows the total number of arcs with length- 
to-width ratio d > 7.5 and d > 10 predicted to be observable 
on the whole sky as a function of <x 8 . Results are shown for 
three different limiting magnitudes in both the I and R bands. 
Arc surveys in X-ray selected cluster samples focused in the 
past on giant arcs, i.e. arcs with length-to-width ratio d > 10 
and R-band magnitudes less than Rn m = 21.5, finding ~ 10 3 
such arcs extrapolated to the whole sky (ILe Fevre et al.ll 1994 : 
Bartelmann et al.lll998l) . 

The numbers given in Figure [10] clearly show how the pre- 
diction falls short of the observation for all values of <x 8 con- 
sidered here, including <x 8 = 0.9. In the latter case, how- 
ever, the difference is only a factor of ~ 2 and can possibly 
be accommod ated including minor contributions due to clus- 
ter galaxies Jlvleneghetti et alJ l2000t iFlores et afl l2000h or a 
slightly different parametrisation for the source-redshift dis- 
tribution. Much progress has bee n made since the first pre- 
diction by iBartelmann et al.l d 1998b . which is also included in 
the Figure. However, if erg is close to the value inferred from 
the WMAP-3 data, all this progress could not alleviate the arc 
statistics problem: with <x 8 ~ 0.75, the predicted number of 
arcs still falls about one order of magnitude below the observed 
number. Should cr% ~ 0.8 persist, there is still a factor ~ 6.5 be- 
tween the prediction and the observation. 



Table 1. Parameters for the two-component Gaussian fit to the 
conditional probability distribution for magnification given a 
threshold d for the length-to-width ratio of simulated gravita- 
tional arcs. 



Parameter 


d = 7.5 


d = 10 


A 


0.84 


0.59 


0-] 


4.1 


4.9 


cr 2 


12.3 


5.8 




11.8 


13.6 




41.3 


38.6 



For completeness, we repeat the calculation of the optical 
depth for arcs with length-to-width ratio d > 10 including the 
scatter in the relation between the mass and the concentration 
of each dark-matter halo in t he population. The p rocedure used 



is the same as described in Fedeli et al. (2007), to which we 



lepr 



refer for details. We just recall that the distribution of concen- 
trations around the nominal value for a given halo mass is well 
fitted by a log-normal di stribution wi t h a rm s of cr c ~ 0.2 ( Jing 
2000). It was shown in iFedeli et alJ (120071 ) that the scatter in 



the mass-concentration relation may increase the total, average 
optical depth by up to 40 - 50%. It was also shown there that 
additional effects like the dependence of the concentration on 
the triaxiality of the dark-matter halos and the distribution of 
projected ellipticities due to the triaxiality itself do not intro- 
duce any additional bias into these results. 

Figure QT| shows the number of arcs with length-to-width 
ratio d > 10 predicted to be observed in a ACDM cosmologi- 
cal model as a function of cr 8 , accounting for cluster mergers, 
with and without the inclusion of the scatter in the halo con- 
centration. To keep the plot readable, we only show results for 
a limiting magnitude of 21.5 in both the I and R bands. The 
scatter in the mass-concentration relation increases the num- 
ber of arcs by ~ 70% for cr 8 = 0.7, and only of ~ 20% for 
cr 8 = 0.9. This decrease is due to the fact that the concentration 
increases with cr 8 , reducing the relative effect of fluctuations 
around the nominal value. On the whole, the scatter in the halo 
concentrations does not help much in improving the agreement 
with observations: models with high <x 8 still fall short by a fac- 
tor ~ 2, and models with low cr 8 are still an order of magnitude 
off. 



6. Summary and Discussion 

We have computed the optical depth, differential optical depth 
and total number of gravitational arcs with length-to-width ra- 
tios d > 7.5 and d > 10 expected in model universes with five 
different normalisations cr 8 . The values for the matter and dark- 
energy density parameters as well as for the Hubble constant 
were taken from the 3-year data release of the WMAP satellite, 
while cr 8 is taken from the set {0.7, 0.75, 0.8, 0.85, 0.9}. 

The cluster population is modelled planting Monte-Carlo 
merger trees for a set of TV = 1,000 dark-matter ha- 
los, uniformly drawn at z — from the mass range 
[l0 14 ,2.5 x lO I5 ]M /i-'. Each halo is assumed to have an 
NFW density profile and elliptically distorted lensing poten- 
tials, with ellipticity e = 0.3. The effect on the strong lensing 
efficiency of the interaction with substructures is also taken into 
account as described in Sect.@] Background sources are prop- 
erly distributed in redshift following to the observed distribu- 
tion of equation |0. 

We converted the average optical depth into an observed 
number of gravitational arcs using an appropriate flux distri- 
bution function for background sources, taking the magnifi- 
cation bias into account. We considered three different limit- 
ing magnitudes in both the I and R bands, including the Ri - 
band limit set by ob servational studies ( Gioia & Luppinoll994 ; 
Luppi no et aDl 1999b . 
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Fig. 10. The number of arcs with length-to-width ratio c/ > 7.5 (left panel) and <f > 10 (right panel) predicted to be observable in 
a ACDM model normalised by erg. Black solid lines refer to I-band, red dot-dashed lines to R-band magnitudes. Results for three 
different limiting magnitudes are shown: 21.5 (botto m pair of curves), 24 (middle pair) and 26 (top pair). The right panel also 
shows the observed number of giant arcs according to Le Fevre et al. ( 1994 ): Gioia & Luppino ( 19941) as a gre en shaded area, to 
be com pared with the lowest red curve. Moreover, the filled black and red points represent the original result of lBartelmann et all 
(1998) (in the I and R bands, respectively), rescaled to the source number counts used in this work. 



Confirming straightforward expectations, we find that the 
total strong-lensing efficiency grows steeply with erg. We find 
that the number of arcs observable in a cosmological model 
with <x g = 0.7 is up to one order of magnitude below the num- 
ber of arcs when cr 8 = 0.9. 

The effect of cluster mergers also depends on cosmology. 
Mergers with relatively small substructures are more likely in 
a low-<x 8 universe, in which individual galaxy clusters are also 
less efficient lenses because of their lower concentration, thus 
causing the effect of mergers to become more pronounced. 
Cluster mergers increase the total optical depth by up to a fac- 
tor of ~ 5 in a model with erg = 0.7 and by a factor of ~ 2 to 3 
in a model with erg = 0.9. 

A particularly strong effect is that the differential optical 
depth at z > 1 for low er 8 can be up to several orders of mag- 
nitude smaller than for high cr 8 . This causes a severe problem 
for explaining the high observed incidence of large arcs in high 
redshift clusters if erg is as low as inferred from the WMAP- 
3 data. Models with low cr 8 also significantly fail to reproduce 
the number of arcs observed in complete, X-ray selected galaxy 
clusters samples. For high <x 8 > 0.9, agreement with the obser- 
vations can be achieved. Thus, we conclude that the arc statis- 
tics problem is unsolved based on the WMAP-3 parameters 
even if a suitable source redshift distribution is included and 
cluster mergers are taken into account. 

We also included the scatter in the relation between 
halo mass and con centration in our calculation, similarly to 
Fedeli et al.l (|2007). This increases the total lensing efficiency 
only slightly, by about 20% for high erg and by ~ 70% for low 
erg. An increment by a factor of ~ 2 may be contributed by 



the inclusion of finer details of the cluster structure, like the 
dark matter halos of sing le member galaxies dMeneghetti et al.l 
2000; iFlores et alJ koOO). The order of magnitude discrepancy 
revealed by models with low erg would however persist. 



We mention that, as we verified, increasing the matter den- 
sity parameter above the fiducial value £2 m ,o = 0.265 used here, 
and keeping the fiat universe assumption, increases the abun- 
dance of observed giant arcs. This is due to the fact that a higher 
£2m,o implies a larger cluster abundance today, and even though 
the evolution of the cluster population is faster due to the minor 
contribution of the cosmological constant, this is not enough to 
counteract the effect in the redshift range important for strong 
lensing. For instance, a WMAP-1 cosmology with erg = 0.9 
and Q m o = 0.3 increases the number of arcs with respect to the 
values reported in Figure [TO] for erg = 0.9, but only of ~ 60%. 
On the other hand, if erg = 0.8 is fixed, the model results can 
be brought in agreement with the observations only with un- 
realistically high values of £2 m o > 0.4. A study of the effect 
on arcs statistics of the combined variation of erg and O m> o, for 
instance following the degeneracy direction given by some cos- 
mological test, is beyond the scope of this work, but is certainly 
interesting for future analysis. 

Some comments regarding the choice of our general as- 
sumption that dark-matter halos have NFW density profiles 
may be in order . Recent publications dNavarro et al. 12004 ; 
Prada et al.l 120061 iMerritt et al.1 120061 iGao et al.l 120071) have 



pointed out that the NFW profile is indeed not the best pos- 
sible representation of dark-matter halo profiles found in high- 
resolution iV-body simulations. This is mainly because the con- 
centrations found depend on the actual radial range used for 
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Fig. 11. The number of arcs with length-to-width ratio d > 10 
and magnitude < 21.5 predicted to be observable in a ACDM 
universe as a function of <x 8 , taking cluster mergers into ac- 
count. Black solid curves refer to I-band, red dot-dashed curves 
to R-band magnitudes. The higher curve of each pair contains 
the scatter in the relation between halo masses and concen- 
trations, while the lower ones are reproduced from Figure [10] 
Again, the green shaded area shows the observational result. 



profile fitting. According to Merritt et al. ( 2006h . the Einastol 
(1965) profile yields better fits, originally proposed to model 
the distribution of stars in the Milky Way. The same authors 
show that the original NFW prescription for relating the con- 
centration and the virial mass of a dark matter halo works bet- 



ter than th e subsequent extensions bv lBullock et al.l ([2001) and 
Eke et alj (1200 lb . even though it must be recalled that they refer 
to the concentrations relative to the scale radius of the Einasto 



profile. However, the lEke et al.l (12.0011) prescription used here is 
sti ll a reasonable fi t to the mass-concentration relation shown 
bv lGao et al. I d2007l) at z ~ 1 , and since there are only a handful 
of efficient clusters at redshifts significantly larger than unity 
even in models with high cr g , we believe this to be a sound 
approximation. 

There are two possible ways of looking at the consequence 
of the results obtained here. Either, the normalisation is high, 
°"8 ^ 0.9, and the low value derived from the WMAP data is 
spurious because of some problem hidden in the data-reduction 
process or in th e Galactic foregrounds (Ide Zotti et al.l 12004: 
Cruz et al. 20061) . Or, and most likely, the actual crg is in fact 
low, as recent weak-lensing measurements seem to confirm 
jFu et al.l 2008), but we fail to properly understand some part 
of galaxy-cluster physics affecting both the number counts and 
the relative strong-lensing properties of galaxy clusters. 

After the submission of this paper, the five-year data release 
of the WMAP satellite was made available. There, the value of 
the matter density parameter slightly increased with respect to 
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Fig. 12. The number of arcs predicted in two examples of cos- 
mological models with early-dark energy, chosen such as to 
agree with CMB, supernovae and large-scale structure data. 
In such models, observational constraints (shown here as the 
green vertical bar) are easily met even with low cr%. Additional 
details on the EDE models presented here can be found in 
Fedeli & Bartelmannl (l2007al) 



the WMAP-3 release, to the best fit Q m , = 0.279 (when CMB 
data are combined with distance measures from la-type super- 
novae and baryon acoustic oscillations). The power-spectrum 
normalization also increased to <x 8 = 0.817, while the error 
bars decreased (0.026 at 68% confidence level for WMAP- 
5). As mentioned above, we expect this slight increase in the 
matter density parameter to be insignificant for our results. On 
the other hand, with crg ~ 0.8, the discrepancy between the 
model predictions and the observations is reduced to a factor 
of > 6. This is less than the discrepancy with WMAP-3 data, 
but still cannot be accommodated by known contributions. For 
instance, different kinds of baryonic physics can increase the 
cluster cross sections of < 100%. In addition, the arc abun- 
dance observed in an X -ray selected cluster sample is only a 
fraction of the total one dFedeli & Bartelmannll2007bl) . This ef- 
fect is not included in the present work, and is likely to worsen 
the agreement with theoretical studies. 

Regarding the abundance of arcs in high-redshift clusters 
(see Figure[7]i, a value of cr 8 ~ 0.8 gives many more distant arcs 
compared to a model with cr g = 0.7-0.75, even though a factor 
of ~ 3 discrepancy compared to high-normalisation models is 
still present. Hence, the high incidence of distant arcs remains 
a problem for a WMAP-5 cosmology, even more than it was 
for a high-normalisation WMAP-1 model. 

In conclusion, the newest WMAP data do not remove the 
discrepancy between theory and observations for the statistics 
of giant arcs, even though it is somewhat reduced with the new 
parameter values. 
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A way out may be possible in presence of a dynam- 
ical dark energy componen t, cf. Figure Q~2] As shown in 



Fedeli & Bartelmann (2007a), the presence of early-dark en- 



ergy (EDE briefly) can play the role of an increased <x 8 on 
non-linear scales because it tends to shift the entire structure - 
formation process to earlier times dBartelmann et al. 2006). 



Consequently, the production of gravitational arcs, and in par- 
ticular the lensing efficiency for high-redshift clusters, are 
significantly increased. Alternatively, non-Gaussian density 
fluctuations may h ave simil arly strong and po s itive effects 



(Gross i et al.l l2007t See also iMathis et alj 12004 ISadeh et al 



120071) . Future studies directed at the recovery of a possible red 
shift evolution of the dark-energy de nsity will be fundamen 



tally important also in th is context dBartelmann et al.l 12003 



Meneghe tti et al. 2005 albl) . Forthcoming analyses of the arc 



statistics problem also require a substantial increase of the ob 
servational data basis, which will be enabled by large-scale op 
tical or infrared surveys in conjunction with fast algorithms 
for automatic arc detectio n (Lenzenetal. 2004; Horesh et al 



20051: ICabanac et alJl2007llSeidel & Bartelmannll2007l) . 
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